Exact asymptotic behavior of magnetic stripe domain arrays 
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The classical problem of magnetic stripe domain behavior in films and plates with uniaxial mag- 
netic anisotropy is addressed. Exact analytical results are derived for the stripe domain widths as 
function of applied perpendicular field, ii, in the regime where the domain period becomes large. 
The stripe period diverges as (H c — ii) _1//2 , where H c is the critical (infinite period) field, an exact 
result confirming a previous conjecture. The magnetization approaches saturation as (H c — H) 1 ^ 2 , 
a behavior which compares excellently with experimental data obtained for a 4 /im thick ferrite gar- 
net film. The exact analytical solution provides a new basis for precise characterization of uniaxial 
magnetic films and plates, illustrated by a simple way to measure the domain wall energy. The 
mathematical approach is applicable for similar analysis of a wide class of systems with competing 
interactions where a stripe domain phase is formed. 



Systems with competing interactions, in particular 
those with short-range attractive and long-range repul- 
sive interactions, commonly develop modulations in the 
order parameter and form domain structures often con- 
sisting of a stripe pattern Realizations are found in 
wide variety of systems, such as magnetic films and plates 
with uniaxial anisotropy [2[, magnetic liquids J3[, type- 
I superconductors in the intermediate state [4|, doped 
Mott insulators [5], quantum Hall structures S, and 
monomolecular amphiphilic ("Langmuir") films [7[. 

The uniaxial magnetic films, where ferrite garnets is 
a classical material studied extensively decades ago for 
use in bubble memory devices, [8j may be regarded as a 
prototype system for stripe domain behavior. Recently, 
the dynamical behavior of the domains in thick garnet 
films showed a vast potential for manipulation of micron- 
sized superparamagnetic beads dispersed in a water layer 
covering the film. By applying magnetic fields with os- 
cillating in- and out-of-plane components, new principles 
for micromachines like colloidal ratchets, size separators, 
micro- tweezers and stirrers, etc. were demonstrated 
Moreover, it has been shown that the magnetic stripe do- 
main structure, when placed adjacent to type-II super- 
conductors, can strongly interact with the vortex matter, 
both in a manipulative way [To| , and as a method to en- 



hance flux pinning in the superconductor [ll|. Thus, one 
sees today considerable renewed interest in the collective 
behavior of magnetic stripe domains. 

On the theoretical side, the treatment of magnetic do- 
mains in plates with perpendicular easy-axis anisotropy 
placed in an external magnetic field is challenging. Even 
solving the magnetostatic problem of one isolated linear 
stripe surrounded by reverse magnetization turned out 
rather complicated analytically, and for a regular array 
of alternating stripes results were so far obtained only by 
numerical calculations [l2|. In this work, based on the 



wall-energy model [2[, i.e., assuming domains separated 
by infinitely thin walls oriented normal to the plate, we 
derive an exact analytical solution for the behavior of a 
periodic array of interacting stripe domains in increasing 
applied field. 

Consider a uniaxial plate of arbitrary thickness, t, 
where magnetic domains form a periodic lattice of par- 
allel stripes with alternating magnetization ±M S , see 
Fig.[TJ In an applied perpendicular field, H, the domains 
magnetized parallel and antiparallel to the field are char- 
acterized by their respective widths and a^, and the 
magnetization of the plate is M = M s (a^ — a^)/a, where 
a = + a± is the period of the stripe lattice. 
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FIG. 1. Cross-section of a plate with magnetic stripe domain 
pattern (top) , and numerical results for the domain widths as 
function of applied field for the case A = 0.05 (bottom). 
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Following the analysis of Kooy and Enz Q , the energy 
density has three contributions; (i) the cost of forming 
domain walls, characterized by the energy cr w per unit 
wall area, (ii) the energy gain of aligning the magne- 
tization with the applied field, —jioHM, and (Hi) the 
self-energy of the domain structure (demagnetization en- 
ergy). The total energy, C7, per unit volume of the plate 
can then be written as 
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Here m = M/M 8 , h = H/M 3 , and A = X/t where 
A = a w /fioM 2 is a characteristic length. The equilib- 
rium magnetization and stripe period at a given applied 
field is given by, dU/dm = dU / da = 0, and expressed by 
the two equations, 
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The particular choice of new variables is motivated by the 
numerical solution of the problem shown graphically in 
the lower panel of Fig. [J only for reference. The striking 
feature is that when the applied field approaches a criti- 
cal value, H c = h c M s where h c < 1, both the period, a, 
of the domain lattice, and the width, a^, of the domains 
magnetized parallel to the field will diverge, whereas the 
reverse domains contract only moderately and terminate 
at a finite width a^ c . At fields above H c , the material re- 
mains single-domain. In the new variables the approach 
towards the critical values corresponds to both y —> 0, 
while the ratio r = x/y = a^/t remains finite [13| . 

Focus of the present analysis is to determine analyt- 
ically the exact behavior as the field approaches the 
critical value. We first derive the relation between h c 
and a^ c . For this, we introduce the auxiliary function 
Pk(z) = Y^Li n~^ kJrl ^ e~ nz = polylog(/c, —z) , and write 

T(x,y) = Im[pi(-ix) - pi(y - ix)]/y. (6) 

For small \z\ one has 

pi (z) = (tt 2 /6) + z(ln z - 1) - (z 2 /4) + (z 3 /72) + . . . , (7) 
which results in the following series expansion, 
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Inserted in Eq. ((2]) it takes the form 
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The critical field is therefore given by 
ft. 



2 r 
— arctan(r~ 1 ) ln(l + r~ 2 ) . 
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(10) 



where r c = a^ c /t is the critical, i. e., the terminal width 
of the minority (anti-parallel to H) domains. 

To find a relation between r c and the material param- 
eter A, a similar treatment is given to Eq. (|3j), using that 
G(x,y) can be expressed as a combination of the real 
parts of both pi(z) and P2(z) with complex arguments 
like those in Eq. (6). For small \z\, one has 



P2(*) = C(3)- 
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where ((n) is the Riemann zeta-function, and Eq. (3) 
becomes 



ln(l + r 2 ) + r 2 ln(l 
270 v 



- 4 )x 4 + . . . = 2ttA . (12) 



The terminal width of the minority domains is therefore 
given by 



ln(l + r 2 c ) + r 2 c ln(l + r~ 2 ) = 2ttA : 



(13) 



and is shown graphically in Fig. [2j The figure also shows 
the dependence h c (A), which follows from Eqs. (10) and 
(13). Both these curves, if replotted as functions of A -1 , 
agree excellently with the numerical solutions presented 




FIG. 2. The critical field, h c = H c /M s , and terminal minority 
domain width, r c , as functions of A = a w /(poMgt). The 
dashed line represents h c (A) = (y/e/7r) exp(— ttA) 
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FIG. 3. Magnetic moment of a ferrite garnet film versus ap- 
plied perpendicular field. The straight line represents the con- 
tribution from the paramagnetic GGG substrate. The inset 
shows a schematic of the two contributions to the moment. 

in Fig. 7 of the Ref. 12i. Note that for any material, 
i. e., given a w and M s , the critical field decreases with 
A, and for A > 0.2 this dependence rapidly approaches 
he = (V^/ 71 ) exp(— 7rA), shown as a dashed line in Fig. 2. 

Consider next the behavior in the vicinity of h = h c . 
Expanding the Eqs. (9) and (12) in Taylor series around 
the critical point one finds to the lowest order that ttx 2 = 
2r c (h c — h). It then follows that the stripe pattern period, 
a/t = r/x, diverges according to 

f-^f (*.-*>-". d4) 

At the same time, the reverse domain approaches its ter- 
minal width as 

a i=r c + : { *:: h L , as) 

t 3 m(l + r c ) 

and the magnetization, m = 1 — 2rt/a, approaches satu- 
ration according to 

m = l-J— (h c -h) 1/2 (16) 

V 7T 

To compare the analytical results with the quantita- 
tive behavior of a typical sample with magnetic stripe 
domains, we prepared a film of bismuth-substituted fer- 
rite garnet, (Y,Lu,Bi)3(FeGa)5 0i2, by liquid phase epi- 
taxial growth on a (111) oriented gadolinium gallium gar- 
net (GGG) substrate. Oxide powders of the constituent 
rare earths, bismuth, iron and gallium, as well as PbO 
and B2O3, were initially melted in a thick- walled plat- 
inum crucible. To ensure homogeneity of the solution a 
stirrer mixed the melt while being kept in the 3-zone re- 
sistive furnace at 1050 °C for 30 minutes. Prior to the 



FIG. 4. Reduced magnetization (M s — M)/M s , experimen- 
tal data and fitted model behavior, Eq.(13), versus applied 
magnetic field. Upper inset: The linear fit to the reduced 
magnetization squared. Lower inset: Magneto-optical image 
showing that the FGF sample displays a parallel stripe do- 
main pattern. 



film growth the melt temperature was reduced to 700 °C. 
The GGG wafer was mounted horizontally in a 3-finger 
platinum holder attached to a shaft rotating by 60 rpm, 
and brought slowly down towards the melt. Finally, the 
substrate was dipped into the melted for 8 minutes re- 
sulting in a macroscopically uniform ferrite garnet film 
(FGF) grown on one side of the substrate. A nearly 
square plate of area A = 21 mm 2 was selected for mea- 
surements. The thickness of the FGF was determined by 
viewing the sample edge-on in a scanning electron micro- 
scope, where a sharp contrast between the film and the 
substrate becomes visible. The ferrite garnet thickness 
was t = (4.0 ± 0.2) fim. 

Shown in Fig. 3 is the result of dc-magnetization mea- 
surements performed using a Quantum Design Magnetic 
Property Measurement System (SQUID magnetometer) 
with the FGF mounted perpendicular to the applied field. 
Above the field of 17 kA/m, the data show linear in- 
crease, which is due to the paramagnetic substrate, in 
combination with the FGF being single domain having a 
constant moment. The fitted straight line intersects the 
vertical axis at a point which determines the saturation 
moment of the FGF sample, ji s = 2.91 10 -6 Am 2 , which 
corresponds to M s =34.6 kA/m. 

Figure 4 shows the reduced magnetic moment of the 
FGF obtained by subtracting the paramagnetic back- 
ground from the raw data. Based on the model result, 
Eq. (16), the reduced magnetization was fitted to the 
predicted asymptotic form (H c — H) 1 / 2 using data over 
a field range below the point where the moment satu- 
rates. The best linear fit of the reduced moment squared 
is seen in the upper inset, from which we find a criti- 



cal field of H c = 13.9 kA/m, and thus h c = 0.40. It 
follows then from Eq. ([TO]) that r c = 0.605, and from 
Eq. CED one finds A = 0.126, and A = 0.126t = 0.504 
microns. The specific wall energy has therefore the value 
a w = 7.58- 10- 4 J/m 2 . 

In previous analyses of the stripe domain problem, see 
e.g., Ref. 12, it was suggested that as the applied field 
approaches h Cl the stripe period diverges with a power 
f) « 0.5. In this work it has been shown that /3 = 1/2 
is an exact result. Consider next what is the field range 
over which the asymptotic behavior is expected to be ob- 
served. There are previous works [HI where experimental 
data were fitted by numerical M — H curves approaching 
saturation seemingly with a finite slope. Furthermore, 
in the classical book Ref. H, the Fig. 2.3 shows M — H 
curves approaching saturation with a finite slope strongly 
depending on the sample thickness. To resolve this ap- 
parent inconsistency, we analyzed the Eqs. (9) and (12) 
up to the next order, i. e., expanding them to (r — r c ) 2 
and keeping the terms ex x 4 . The analysis shows that 
Eqs. (14)-(16) provide a good description as long as 
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(h c -h) < min{r c , 1} . 



(17) 



Thus, for thick plates, t ^> A, one has a very small 27rA 
and r c ~ [27rA/ ln(l/27rA)] 1 / 2 , which is much less than 
unity. Thus, it follows that the asymptotic behavior (16) 
will be observed only very close to h c . In practice, the 
field interval may be beyond experimental resolution, and 
the slope of the magnetization curve near h = h c appears 
finite. For thinner plates and films, r c rapidly increases, 
see Fig. 2, and the inequality (17) becomes much weaker, 
and the range where one should observe the critical be- 
havior Eq. (16) will be sizable, as demonstrated in the 
present experiments. 

Finally, note the presence of a small shoulder in the 
reduced magnetization data in Fig. 4 seen just above H c . 
Here the behavior deviates significantly from the model 
prediction, and is caused by topological fluctuations in 
the domain pattern. The role of material defects caus- 
ing wall pinning becomes significant, and the stripes lose 
their alignment, a scenario which is readily seen visually 
by following the behavior using magneto-optical imag- 
ing 0. 

In summary, we have presented an analytical asymp- 
totic solution to the problem of modeling the behavior of 
an infinite array of parallel alternating magnetic stripe 
domains subjected to a transverse field. The mathemat- 
ical approach used in this work can be applied to derive 
exact results also for other systems [H, [Hj where stripe 
domain phases are formed and described by a conflgu- 
rational energy term similar to the one treated in the 
present case. 
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